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Abstract 

The symmetric and traceless part of the matrix element ST{N\qD^xDyq\N) can be determined from 
the second moment of the twist-3 parton distribution function e{x). Recently, novel experimental data 
on e{x) have become available, which enables us to evaluate the magnitude of the above matrix element 
with considerably reduced systematic uncertainties. Based on the new experimental data, we show that 
ST{N\qD^Dvq\N) is likely to be at least an order of magnitude smaller than what previous model-based 
estimates have so far suggested. We discuss the consequences of this observation for the analysis of deep 
inelastic scattering and QCD sum rules studies at finite density for the vector meson and the nucleon, in 
which this matrix element is being used as an input parameter. 
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I. INTRODUCTION 


The traceless and symmetric component of the matrix element of the operator qD^Dyq between 
a one-nucleon state, 


ST{N\qD^Dyq\N) =\sT[ {N\uD^Dyu\N) -h {N\dD^Dyd\N)], 

1 


= -e2[PiiPv-^glivM] 


( 1 ) 


TV u 


is known to play a role in the analysis of dera inelastic scattering (DIS) yj] and in applications of 
the QCD sum rule method to finite density yllj]. In Eq.O, stands for the four-momentum of 
the nucleon and is the nucleon mass. The value of € 2 , however, has so far only been estimated 
via simple models or certain assumptions on the proportionality between different matrix elements 
[PD. While such estimates may be fine for obtaining a first qualitative idea on the magnitude of € 2 , 
it is far from clear whether they are quantitatively reliable. 

The situation now changed with the availability of new experimental data [|4l], which, as we 
will see, strongly constrain the value of 62 - This is possible due to the fact that €2 is related to the 
second moment of the twist-3 distribution function e(x) as shown below: 


€ 2=1 dxx^e{x) 

Jo 

= - [ dxx^\e‘‘{x)+e'^{x)+e“{x) + e‘^{x)]. 
2 70 

Here, the various flavor components (x) are defined as [qI] 


( 2 ) 


e^x) = 


1 


2Mn 




(3) 


with [0, Xn] being the gauge link for making the above expression gauge invariant and n a null 
vector with mass dimension —1. Through the analysis of experimental data on the beam-spin 
^mmetry of di-hadron semi-inclusive DIS obtained at the CLAS experiment at Jefferson Lab 
jdll , it has become possible to extract a small number of data points for e^{x), which is defined as 
follows: 


' W = ^ [e"{x) - e"(x)] - ^ [e‘^{x) - e^(x)]. 


Making use of some reasonable assumptions on the flavor structure of e{x) and on its behavior in 
those X regions, where no data points are available, will allow us to get an estimate of €2- 

As a result, we find that even though the experimental uncertainties are still rather large, the data 
can constrain the magnitude of €2 to be at least an order of magnitude smaller than values obtained 
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from the previous simple estimates yy. This means that the matrix element ST{N\qD^Dvq\N) 
has been largely overestimated in the past DIS or QCD sum rule analyses at finite density. 

Indeed, applying the novel estimate of 62 to the operator produet expansion (OPE) of the elee- 
tromagnetic current, which contains information on the spin-averaged structure functions F2 and 
Fi, it is found that the above matrix element only gives a contribution of 3 % or less compared to 
the experimentally extracted values of the twist -4 effects in the second moments of these structure 
functions and can, therefore, be ignored at the presently available level of precision, which is in 
contrast to the conclusions of earlier studies. Furthermore, examining the OPE of the vector cur¬ 
rent correlator in nuclear matter, coupling to the p, CO, and (j) mesons, and, separately, the nuclear 
correlator in nuclear matter, we similarly find the relevant contributions to be small. 

The paper is organized as follows. After explaining how to extract the value of €2 from the 
experimental data in Sec. HI we study its consequences in Sec. IIIl which includes a discussion of 
both the OPE needed for analyzing DIS data and for the sum rule analyses of vector mesons and 
the nucleon at finite density. The paper is summarized and concluded in Sec. UVl 


II. ESTIMATION OF ^2 
A. Earlier simple estimates 

Before discussing the estimation of 62 based on the newly available experimental data, we here 
for illustration and later comparison briefly describe two simple methods, that have so far been 
used to compute €2- 


1. Method 1 


Here, we will follow yj] to estimate the needed matrix element. 


{N\qD^Dyq\N)c^-P^P^{N\qq\N) 

1 _ ( 5 ) 

= -^PliPv{N\qq\N). 

P^ in the first line represents the average momentum of the quark q in the nucleon, while the 
second line follows from the assumption that P^ is about 1/6 of the nucleon momentum p^, as 
half of the nucleon momentum is carried by the gluons and the rest is divided evenly among the 
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three valenee quarks. Transforming this result into a traeeless form, we get 


ST{N\qDiiDyq\N) ~ -^{N\qq\N){p^Pv-\giivMf,), 


( 6 ) 


whieh hence means: 


^2 ^ ^(iV|^^|A^). 


(V) 


Expressing this through the quark mass + m^) and the kN sigma term cr^jiv = 

m„(A^|^^|A^), we get 


62 ^ 


1 OjiM 

36 m„ 


( 8 ) 


Using CT^riv = 45MeV [a] and niq = 3.5MeV DTI], we finally obtain 


62 — 0.36. 


(9) 


2. Method 2 


Here, we briefly recapitulate the c 
Combining Eqs.(4.43) and (4.49) of \2\ 


iscussion of yfl to estimate the needed matrix elements. 
, we obtain 


ST{N\qDaDpq\N) = ^^(^{N\qDQDoq\N) - ^{N\qgO ■ Gq\N)^ (^PaPh - , (10) 

where the second term follows from the fact that the two covariant derivatives do not commute. 
The above equation then immediately gives 


4 1 


3m2 


1 


^2 = -tttt (A^|^E>o^o^|A^) - ^{N\m(y ■ Gq\N) . 
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( 11 ) 


Employing a bag model estimate, it was furthermore shown in yO that first term in the bracket on 


the right-hand side of the above equation is much smaller than the second one and can, therefore, 
be ignored. 

Next, using a parametrization proposed in [|2|l 


(A^|^gC7-G^|A^) =ml{N\qq\N), 


( 12 ) 


we are led to 


1 , , 


1 ^0 

6mI niq ■ 


( 13 ) 
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mp is believed to be about ~ 0.8 ±0.2GeV^ in vacuum [Isl], which was assumed in yD to hold 


also for the one-nucleon state. This assumptions leads to 


£2^1.95, (14) 

which is about 5 times larger than Eq.® of the previous subsection. 

As we will see later, the estimate in Eq. (fT4l) turns out to be 2 orders of magnitude larger than 
our updated value based on experimental constraints. This means that the ansatz of Eq. lfT^ with 
Mq ~ 0.8 ±0.2GeV^ is most likely an overestimation of the actual matrix element. We, therefore, 
caution the practitioners of QCD sum rules at finite density to be careful when making use of this 
parametrization. 


B. Evaluation based on experimental data 


As already mentioned in the Introduction, the only information that we presently have from 
experiment, is the magnitude of {x) [given in Eq.®] at a few values of x. This, obviously, does 
not suffice to determine 62 completely and we, thus, have to introduce a number of assumptions 
on the relations of e“{x), e^{x) and their sea-quark counterparts. Eor this, we will have to rely 
partially on model calculations of e{x). i^ecifically, these models are the bag model (BM) [SI], 
the chiral quark soliton model (;tQSM) [| 9 l], and the spectator model (SM) 1I10I1 . To get an idea 


of the systematic uncertainties of these assumptions and models, we will test several versions of 
them and study their effects on 62 - Eor illustration we give the 62 values and their respective flavor 
decompositions obtained from the BM, ;t^QSM, and SM in Table H 

Throughout our whole study, we will assume the sea-quark effect on 62 to be flavor symmetric. 
Namely, we will set 

ri _ /■! _ _ 

(15) 


[ dxx^e'*{x) = [ dxx^e^{x) = e^- 

h Jo 


The violation of this flavor symmetry can be studied by both the BM and ;t^QSM, which show that 
it is only a very small effect (see Table U), which can be ignored here. 

Next, we will have to fix the relative strength of the u- and J-quark contributions to the second 
moment of (x). One could naturally expect that the two are proportional to the number of 
respective valence quarks; hence. 


^2 = J dxx^e‘^{x) = ^ J dxx^e“{x) = ^^ 2 : 


( 16 ) 
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TABLE I. Second momentsjof e(x) for the various quark components, extracted from the bag model [SI], the 

H. 


chiral quark soliton model 


and the spectator model 



BM 

ZQSM 

SM 

el X 10^ 

8.1 

5.3 

10.7 

X 10^ 

4.1 

3.1 

-2.2 

el X 10^ 

0.2 

0.4 

- 

el X 102 

0.3 

0.3 

- 

C2 X 10^ 

6.3 

4.6 

4.3 


which is satisfied with good accuracy by both the BM and ;^QSM. We will call this assumption 
“Ansatz 1” in the following. 

The SM, however, seems to suggest a somewhat different picture, in which e^(x) shows an 
oscillating behavior and its second moment, therefore, even becomes a negative, but rather small 
number (see Fig. 9 of 111 ill and TableH])- We will, therefore, set as a second assumption 

/' 


dxx^e‘^{x) = 0, 


(17) 


which we call “Ansatz 2”. 

For Ansatz 1, we can rewrite the second moment of the experimentally measured function (v) 
as 

[ dxx^e^{x) = ^€2 - \e 2 , (18) 

Jo to 3 

and for Ansatz 2, 

[ dxx^e^{x) = 3^2 - \e 2 . (19) 

Jo 9 3 

As a last point, we need to fix the the ratio between the second moment of the u and u quarks, 
for which we can obtain some guidance from the BM and the ;kQSM. For the BM, the e"{x) is 

small and has some sizable strength only around x = 0. Its second moment is, hence, negligible. 

On the other hand, the e“(x) for the ;tQSM is larger, with having the size of almost 10 % of the 
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valence value of e\. We, thus, define 


dxx^e“{x) = rf I dxx^e'^{x), 


2 Ml 


( 20 ) 


and choose for the parameter r\ a range of 0 < r/ <0.10. 

With this definition, we can finally relate the second moment of the experimentally measured 
function {x). 


dxx^e^ {x) = 62 , 


( 21 ) 


with ^ 2 ,^ 2 , ^2 and e^- For Ansatz 1, we get 


18 


^2 = 


61 = 


6^=60 = 


7-6t] 

9 

7-6t] 

18?7 




'2 5 


( 22 ) 


7-6t] ^ 

The results for Ansatz 2 can meanwhile be given as follows: 


^2 = 


4-3?7 


4, 


ei = 0, 


(23) 


el = ei = 


9t? V 


4-3t] 

Next, we briefly discuss how to obtain from the experimental data of In the region 
of the three data points provided by experiment, we employ the most simple rectangular rule to 
approximate the integral, as indicated in Fig. [B For applying this rule, the horizontal rectangle 
sizes are fixed as follows: 1) The boundary of any two neighboring rectangles is determined to 
be in the center of the corresponding data points. 2) The lowest and highest rectangles are set 
to be symmetric with respect to their data points. For the region below the lowest rectangle, 
we make use of the fact that (x) at small x is known to have a behavior with a < 

1 [SI]. Note that the leading Pomeron contribution (a = 1) vanishes here because quarks and 
antiquarks contribute with opposite signs to ^ (v). However, as we have no additional information 
on the value of a, we will nevertheless set a = 1 as an upper limit, fix the coefficient A from 
the lowest data point and compute the respective contribution to the second moment analytically. 
Furthermore, the experimental points appear to quickly approach 0 above the highest data point 
at X = 0.356. We, thus, assume it to be 0 for x values above the respective highest rectangle. 
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FIG. 1. (Color Online) The three experimentally available data points in both the WW scenario (red points) 
and leading scenario (blue points), shown together with the rectangles used for the numerical integration, 
as described in the text. In addition to the data points, the extrapolation used for {x) at small x, for which 
we have assumed a 1 /x^-type form, is illustrated below x ~ 0.13 as red and blue lines, respectively. 


All this is pictorially illustrated in FigUl in which the shaded areas depict the way the numerical 
integration is performed. 


We should mention here that computing the second moment of e{x) with only three available 
data points is in principle an ill-defined task. However, all of the models describing the function 
e{x) have their dominant strength below 0.5, as is shown in Figl2l Moreover, for the small-x 
behavior of e(x), we are using a conservative upper limit, which is based on general considerations. 
Therefore, even though we have to rely only on three data points, the values of e{x) at these 
points can be expected to determine the order of magnitude of its second moment, which is all the 
precision needed for the present work. 


The experimental results of M, in fact, contain the data of two different analyses. The first one 
is based on the Wandzura-Wilzcek approximation (and is, therefore, called the “WW scenario”), 
while the second one includes terms that go beyond the Wandzura-Wilzcek approximation and is 
called the “leading scenario”. In the following, we will in the following consider the results of 
both scenarios to get a rough estimate of the systematic uncertainties involved. The numerical 
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FIG. 2. (Color Online) The function e{x) for the three models (BM, ;^QSM and SM) considered in this 
work. The numerical data needed for this plot have been extracted from Fig. 9 of | 


im 


integration described above and depicted in Fig [T| then gives 


"1 I (2.35 ±0.32) X 10-2 (WW scenario), 

I dxx^e^{x) = 

® ’ (—0.97 ±0.32) X 10-2 (leading scenario) 


(24) 


All components are now in place for estimating 62 - First of all, using Eq.© we can combine 
the various flavor contributions given in Eqs. (l22l) and (l23l) as 

1 


^2 — 2 (^2 + ^2 + ^2 + ^2) 


(25) 


Applying then the results of of both the WW and leading scenario, we get altogether four 
values of 62 with ranges determined by the variation of 7 ]. These ranges are shown in Table 
mi together with the corresponding flavor decompositions. It can be seen in this table that the 
largest uncertainty of €2 is related to the discrepancy of the two scenarios used for analyzing the 
experimental data. Comparing the numbers of Table HIl with those of Table H it is observed that 
the model values have the same order of magnitude, but are generally somewhat larger than the 
experimentally extracted ones. Eurthermore, it is noted that model values are mostly positive, 
which means that they favor the WW over the leading scenario (see also the recent discussion in 

Imx 

With the four ranges obtained above, we can now give our final estimate for e 2 - Taking the 
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TABLE II. Second moments of e{x) for the various quark components, extracted from the experimental 
values of 1^. The ranges in each entry of the table are obtained by employing the central values of Eq. (l2^ 
and varying the pai'ameter t] of Eq.(l20l) between 0 and 0.10. 


WW scenario leading scenario 



Ansatz 1 

Ansatz 2 

Ansatz 1 

Ansatz 2 

c" X 10^ 

6.0 ~ 6.6 

5.3 ~ 5.7 

-2.5 ~ -2.7 

-2.2 ~ -2.4 

cf X 10^ 

3.0 ~ 3.3 

0.0 

-1.3-1.4 

0.0 

= X10^ 

0.0 ~ 0.7 

0.0 ~ 0.6 

0.0 ~ -0.3 

0.0 ~ -0.2 

62 X 10^ 

4.5 ~ 5.6 

2.7 ~ 3.4 

-1.9 ~ -2.3 

-1.1 ~ -1.4 


smallest and largest value in the bottom line of Table HI] and, furthermore, allowing for the possi¬ 
bility that ^2 can vary in the range determined by the statistieal errors given in Eq. (l24l) . our result 
reads 

e2 = -0.030 ~ 0.064, (26) 

whieh should be eompared to the numbers 0.36 of Method 1 and 1.95 of Method 2. This eompari- 
son elearly shows that both Methods 1 and 2 have overestimated the relevant matrix element by at 
least an order of magnitude. Therefore, even though the experimental uncertainties of 62 are still 
large, we can be quite certain that it must be much smaller than previously expected. 


III. DISCUSSION OF EFFECTS ON DIS AND QCD SUM RULE ANALYSES 
A. Contribution to deep inelastic electron scattering 

The relevant OPE can be obtained by using the electromagnetic current j^{x) = q{x)QY^q{x) 
in the correlator 

n^y{C 0 ,q) = i I j/e'^"(iV|T[j^(x)yv(0)]|iV). (27) 
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Here, Q is the charge operator, and \N) stands for t 
including the target mass corrections, appear as dl] 


le one-nucleon state. The twist-4 contributions. 


n 


twist 4 1 

jUV ' 


CO,q) = 




1 


13 


d^v ( + g ^ —-^C + E ] + e 


4^ 2^ 


9 Q 


1111 
(p.,)2V4'*+4®“2'^+2^ 


(28) 


where we have defined 


and 


— Sp.v 


^2 ’ 


dpLv — 


P^Pv 

{p-qf 


q^ + {p^iqv+Pvq^) - g^iv 


p-q 


(29) 

(30) 





2 {p-q)' 


(31) 

(32) 


In the above definitions, stands for the four-momentum of the nucleon, with Ma? being the 

nucleon mass; hence, The parameters A-E are related to the operators shown below: 


= gSTq Gp^] YaQ^q = g^ST{qYat‘'Q^q) Y^^qY^t^q), 

q 

Bali = gSTq{iDa,Gp^}Y^Y5Q^q, 

Cap =mqSTqDaDpQ^q, 

B>ap = g^d'q[Da,G^p]Y^lQ^q, 

Eap =g^ST{qCY5YaQq){qt''Y5YpQq)- 

The expectation values of these operators are parametrized as 

{N\Aap\N) = [paPp - ^Mlg^p^A, 

WapW) = (paPp - ^Mlgap^B, 

{N\Cap \N) = (paPp - ^M^gap^, (34) 

{N\Dap\N) = [paPp - ^Mlgap^D, 

{N\Eap\N) = (paPp - ^Mlgap^E. 
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In the deep inelastie limit and x = finite, one ean negleet the target mass eorreetions, and 

the matrix elements eontribute to the seeond moments of the transverse and longitudinal strueture 
function as follows: 


F^='^dx = 0.005 ± 0.004 GeV^(proton) 




13 


(35) 


C), 


C Fl=^dx = 0.035 ± 0.004 GeV^ (proton) 

Jo 


= 2M,(Ia-Fb-Ic). 


(36) 


Here, the values in the first lines are obtained from experimental constraints as extracted in [J]. 
The overall factor of 2M]\f appearing in the second lines comes from the covariant normalization 
factor of the nucleon used in extracting the above numbers. Using the notation introduced before, 
the value of C appearing in Eqs.(l35]) and (l36l) can be written as 

(37) 




Using the values given in Table|II]and the PDG averages for the quark masses at 1 GeV [Tl], we find 
that the range of C values are given as in Table Hill As one can see from the table, C contributes 
less than 4% and 1%, to the extracted numbers of the second moments for the transverse and 
longitudinal parts, respectively, for all cases. Hence its contribution can be safely neglected. It 
should be emphasized here that this conclusion differs from yl, in which an old estimate in line 
with Eq.® was used, and where the contribution of C to the second moment of the transverse 
structure function was found to be sizable. The novel findings of the present work have, therefore, 
somewhat changed this situation. 


B. Contribution to the p and ft) meson sum rule 

The relevant OPE to the vector meson sum rule can be obtained by using the vector meson 
current j^{x) = q{x)T:Y^q{x), where T is the isospin operator, in the correlator 

n^y{C0,q) = i I j/e'^"(T[7^(x)7v(0)])p. (38) 

Here, ()p stands for the expectation value with respect to the ground state of nuclear matter at 
r = 0. In the linear density approximation, the matrix element can be obtained as the nucleon 
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TABLE III. Values of C that contribute to the second moments of the structure functions in units of GeV^. 
The ranges in each entry of the table are obtained by employing the central values of Eq. d^ and varying 
the parameter rj of Eq. (l20l) between 0 and 0.10. The second and third lines use factors as they appear in 
Eq.d^ and Eq.d^. respectively. 



WW scenario 

leading 

scenario 


Ansatz 1 

Ansatz 2 

Ansatz 1 

Ansatz 2 

2MnC X 10^ 

-9.9 ~ -12 

-6.8 ~ -8.5 

4.1 ~ 5.0 

2.8 ~ 3.6 

2Mn^C X 10^ 

16~20 

11 ~ 14 

-6.7-8.1 

-4.6 ~ -5.9 

2Mn^Cx10^ 

22~27 

15 ~ 19 

-9.2 - 11 

-6.3 ~ -8.1 


matrix element times the density p. Hence, the twist-4 operators contribute to the OPE as in 
Eq. (l2^ multiplied by p, the only difference being that the operators in Eq. (l3^ contain the isospin 
matrix instead of the charge operator Q. In the sum rule analysis, the vector meson is taken to be 
at rest ^ = 0, in which case there is only one invariant tensor, and the Borel transform is taken with 
respect to — —)■ oo. The sum rule including the twist-4 operators, except the operator of interest 


here, namely , are given in [[13 


m. 


While the matrix elements are defined with the charge operator Q replaced by the isospin oper¬ 
ator T, the relative contribution of the operator 0^^ relative to the other twist-4 operators remains 
small as in the previous subsection, and its contribution to the vector meson sum rule can be 
safely neglected. The situation could, however, be different for the (j) meson as the small mu.md 
is replaced by the larger strange quark mass m^. This case will, thus, be considered next in more 
detail. 


C. The OPE for the (j) meson channel in nuclear matter 


The 0 meson can be described by the interpolating field j^{x) = s{x)y^s{x), which is substi¬ 
tuted into the two-point function of Eq.dSS]). Up to dimension six and twist-2 terms, the OPE 


for this correlator has already been given in earlier works yj 


16h . Here, we are interested in 
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the operators of dimension six and twist-4 that contain strange quark fields. Within the linear 
density approximation, the Wilson coefficients of such operators are obtained in analogy to the 
previous subsections. The final result can then be given as Eq.(l28l) multiplied by p, the operators 
corresponding to the parameters A-E now being 


gST 5 JaS = g^sr[sjafs) 


Bap = g^Ts{iDa,Gii^}Y^r5S, 
Cap =msSTsDaDps, 

Dap = g^Ts[Da,G^p]r^s, 

Eap = g^^'CisCr5ras){sCY5rps). 


(39) 


It is noted that, as before, D does not appear in the final result and we, therefore, do not need to 
be concerned with the corresponding operator any longer. An estimate of the one-nucleon matrix 
elements of the other operators is given below. 

1. Numerical estimates of the twist-4 matrix elements 

We here wish to evaluate the one-nucleon matrix elements of Eq. (l39l) . This can, however, not 
be done as reliably as for the u- or J-quark case, because no experimental information from deep 
inelastic scattering is available. Nevertheless, one can still try to get an estimate by making certain 
assumptions on the behavior of the matrix elements. Specifically, we will make repeated use of 
the ansatz 


{N\srOs\N) ~ {N\urOu\N)^, 


(40) 


for relating the i’-quark operators with a general operator insertion O to their w-quark counterparts. 
Here, the parameters Aj and A" are moments of parton distributions of the nucleon, with quark 
flavor s and u. They are defined as 



(41) 


( 42 ) 
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Using the recent estimation of the parton distribution functions in lll7l] . the moments can be eval¬ 
uated as 


A'; = 0.808 ±0.069, 

Aj =0.0443 ±0.0102. 


(43) 

(44) 


ST(N\g^(-syal‘‘s)Y.,(qyiit‘‘qm 

Using Eq. (l40l) . we rewrite the matrix element as follows: 

ST{N\g\syat‘^s) Y,{qy^fq) \N) ^ ST {N\g\uyafu) '£iqypt‘^q) |iV)§. (45) 

q q 

Ignoring the strange quark contribution to the sum on the right-hand side of this equation, it 
is seen that that it can be related to the parameter K^, discussed in yfl: 


1 


1 


.2^ 

All- 


ST{N\g^{syat‘'s) ^(^r/3^"9') \N) ^ {paPp - -Mlgap)K-, 

Therefore, the parameter A can be given as 

1 2^1 
A = _ 

2 Mn “Af 

ST{N\g^iDa,Gp^}y^y5s\N) 

The next matrix element can be treated in a similar way: 

Sr{N\g4iDa,Gii^}y^y5s\N)c^Sr{N\u{iDa,Gii^}y^y5u\N)^. 

- — 


(46) 


(47) 


(48) 


Using Eq.(19) of yfl, which is rewritten as Eqs.(Cl-3) in [Il8h . it can be related to the pa¬ 
rameter Ki: 


1 


Sr{N\gs{iDa,Gp^}y^y5s\N) - :^{PaPp - 


1 


Comparing the above result with Eq. (l34l) . we can express B as 

1 

B = ^KS^. 


'A“ 


2Mn ^“Af 


(49) 


(50) 


msST{N\sDaDps\N) 

To study this matrix element, we will make use of the knowledge of the previous sections, 
which dealt with {N\qDaDpq\N). Using again Eq.(l40l). we get 


m,ST{N\WaDps\N) - m,ST{N\uDaDpu\N)^ 


= -inse2 


2Mn 


2 „ \ ^ 1 


PaPp - 


(51) 


Af 
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which means that C can be given as 


C = -mse2 


A\ 
Af 


(52) 


ST {N\g^{st‘^ysyas) (st^’ysyps) \N) 

The last matrix element to be determined is the one eontaining four strange quarks. We use 
the same strategy as above and rewrite it as 

ST{N\g^{sty5yas){st‘‘y5yps)\N)c^Sr{N\g^{sty5yas)iut‘‘y5ypu)\N)^ 


~ ST{N\g^{ufy5yau){ut"y5ypu)\N) 


AjN2 

A" 


(53) 


2Mn 


1 \ 1 
PaPp - ^M^gap {Ku - 


A" 


In the last line, we have used Eqs.(Cl) and (C3) of rtl8h . We can, hence, express the param¬ 
eter E as 


E = 


2M^y “ 2 


(54) 


For extraeting the values of A, B, C, and E from the formulas of the this subseetion, we need 
the values of the parameters K^, K^, Ki, and These were estimated in yfl from lepton-hadron 
deep inelastie seattering measurement results and by assuming ratio of K^, and K^d to lie in a 
reasonable range. We will use the findings of yO together with the PDG average of the strange 
quark mass at 1 GeV DTI] to give an estimate of our parameters. The eorresponding result is given 
in Table HVl 


2. Consequences for the (j) meson sum rules at rest 

To get an idea of the size of the eontributions of the twist-4 terms, we will here evaluate their 
eontribution to the OPE of the ^ meson at rest. A more detailed study on the OPE for the general 
finite three-momentum case will be done elsewhere. 

Eet us simplify Eq. (l28]) by setting p = (Myy,0,0,0), whieh means that we are working in the 
nuclear matter rest frame. Furthermore, we set q = 0, take the traee over the Eorentz indices and 
define 

n(^^) = -:^nJ(^o,^ = 0), (55) 
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TABLE IV. Values of the parameters A, B, C, and E in units of MeV, for six different data sets of KI, K^, Kf, 
and given in yl]. For C and Coid, we use Eq. (l5^ with the central value of our new estimate of Eq.d 
^2 — 0.01, and the old value of Eq. (fT4]) . respectively. 




A 

B 

C 

Cold 

E 

Kl-- 

= 0(5 

5.93 

-6.95 

-0.12 

-13.69 

-1.56 

K-- 

= C(/5 + l)//5 

3.21 

-8.76 

-0.12 

-13.69 

0.22 

Kt = 

= kL 

1.93 

-9.61 

-0.12 

-13.69 

1.07 

Kl-- 

1 

II 

-5.29 

-14.43 

-0.12 

-13.69 

5.91 

Kl- 

= -KUI3 + 1)/I3 

-6.57 

-15.27 

-0.12 

-13.69 

6.76 

K-- 

= 

-9.29 

-17.08 

-0.12 

-13.69 

8.54 


which gives for the dimension-six, twist-4 term, 

ndlm.6,,w«4(^2) ^ 4MjpL ( -Ia-Ib-Ic-\e 

% 

Considering the expression of X together with the numbers of Table |IVl it is seen that B gives 
the dominant eontribution for all data sets and that C with our new estimate is negligibly small. 
It is, however, noted that with the old estimate of Eq. (fT4l) . whieh would mean using Com instead 
of C, the C-term would, due to its large prefaetor 7/4, aetually have provided by far the largest 
eontribution and would have inereased the overall value of X at least by a faetor of about 4. If the 
findings of this paper are eorreet, this will help to largely reduee the twist-4 eontributions to this 
sum rule. 

Making use now of our new estimate of C and of the other parameters given in Table UVl the 
variable X is found to take values that lie roughly between 4 and 9 MeV. After Borel-transforming 
the result of Eq.(l5^. we obtain 



■p-jdim.6,twist4 

^^OPE 




( 57 ) 
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Let us compare this finding with the other terms appearing in the OPE at the same order. These 
are given for instance in Eq.(12) of [Il6n . from which we list below the linear density coefficients 
of the dimension-six term: 



(58) 


(59) 


Numerically, these give 


1.71 X lO^MeV^-S.SOx lO^MeV^ + 2.01 x lO^MeV^ 
-4.19 X lO^MeV^-2.80 X 10^MeV^-7.61 x 10^MeV^ 


(60) 


(61) 


for which we have used K)v = 1 and Osn = tns{N\ss\N) = 50MeV, which is rather on the larger 
side considering the recent lattice QCD calculations of this quantity. Eor the other parameters, we 
used the values given in ll^ . It is clear from the above numbers that the four-quark condensate 
term is dominant. All this now has to be compared to Eq.dSTI). without the density and Borel mass 
factors. This leads to 


-2XM^~ -1.13 X 10^MeV^ 


(62) 


where we have set X = 6.40MeV, which lies at the center of the expected value range of X. 


Even though our new estimate of the C operator has helped to reduce the magnitude of the twist-4 


contribution, the result of Eq. (l62l) is still rather large, namely, of the same order as four-quark 
condensate term and more than 16 times larger than the twist-2 term. 

D, Consequence on the nucleon sum rule and nuclear symmetry energy 


The sum rule for the proton can be obtained by looking at the correlation function of the proton 
interpolating fields. 




= ns{q^,qu)+nq{q^,qu) ^ + Tlu{q^,qu) 


(63) 


where q is the external momentum, u is the medium four-velocity and Y\s{q^,qu), Y\q{q^,qu), and 


Ylu{q^Tqu) are the three invariants according to Eorentz, parity and time reversal symmetries [Il4l] . 


The Ioffe current [l20n will be used as a proton interpolating field: 


ri(x) = £ahclu^(x)CY^Uh(x)]Y5Y^dc(x). 


(64) 
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1. OPE for the proton 


The OPE for proton correlator of Eq. (l6^ is given, for instance, i n lllSh . Including, furthermore. 


the Us corrections 1121 
follows: 

tOPE/ 2 


2311 and the dimension-eight contributions 


Imi 


it can be summarized as 


^(-^^)ln(-^VM^) + {dd)pj + ^^{d{iDQiDQ]d)p^i, (65) 






+ 


1 


1 


6 \ q 

+ (qu) 


{uu){gsUoQu)pj 




{{u^u)pj + {d^d)pj), 

( 66 ) 


n^^\q\qu)^- 


(67) 


Here, /i is the renormalization scale which will be matched with the Borel mass M after the Borel 


transformation. Note that the operator of interest appears as the 


ast term in the scalar self-energy. 


The density dependencies of all the operators are given in lllSh except for that of dimension 
eight, for which we use the approximation 


{[qq]){[gsq<y-Qq])p,i = {[qq]h^c{[gsq(y ■Gq])^^c 

+ {{[qq]) vac([gj^Ct • Qq\Q)p -|- ([^^]o) p(\_gsq^ ■ Gq\)\a.c) 




n+{mq) 


p, 

( 68 ) 


where ■ ^<2'])vac = (0-8)([<?^])vac and and ” stand for u and d quark flavors, respec¬ 

tively. {[gsqo ■ Gq]Q)p is chosen to be 3 GeV^ as in 118.1251]. The detailed definition and description 
of the ratio TZ±{mq) can be found in Sec. Ill B of lll8n . 


2. Result for the nucleon sum rule 


The analysis for the nucleon and the symmetry energy follows that given in [Il8h . with / = —0.3 


Here, / parametrizes the degree of factorization of four quark operator in medium {qq)j, (1 
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FIG. 3. (Color Online) Borel window for ;B[n^(< 7 Q, |^|)] (a) with OPE up to dimension-six condensates 
and (b) dimension-eight condensates together with corrections. The solid (dashed) lines show the ratios 
between the contribution of the highest dimensional operators (the continuum) and that of the total OPE. 
The horizontal dotted lines correspond to the ratios being 0.75 (red) and 0.5 (blue) in (a) and 0.05 (red) and 
0.5 (blue) in (b), respectively. 


f){w)lac + /(w)n- The formalism and other parameters are based on the nueleon sum rule at 


> oo but fixed \q\ 


[Q. 


First, we show the Borel window before (a) and after (b) ineluding the eorreetions and 
dimension-eight eondensates in Fig. [3l The left panel of Fig (3) is adapted from (IlSll . There, the 
upper bound for the Borel mass was obtained by requiring the ratio between the eontribution of the 
highest dimensional operators and that of the total OPE to be less than 50%, and the lower bound 
from the eondition that the ratio between the eontinuum eontribution and the total OPE should 
be less than 75%, whieh eonfines the Borel window to 1.0 GeV^ < m 2 < 1.2 GeV2. The added 
eorreetions lead to a eonsiderable improvement, as ean be inferred by eomparing plots (a) and (b). 
We ean, therefore, now require above the ratios to be less than 50 % and 5 %, respeetively, whieh 
eonfines the Borel window to 0.8 GeV^ < M^ <1.0 GeV^. 

Next, we study the sum rule for the sealar (M*), veetor self-energy (Ey), and the ehange in 
the pole mass [Eq{I = 0)] for symmetrie nuelear matter, whieh ean be obtained from eombining 
the three invariant funetions appearing in Eq. d^ as diseussed in lllSn . As ean be seen in Eig. |4] 
(a), the ratio £’g(/=o)/^A/^ varies from 0.90 to 0.84 and M^/Myy from 0.55 to 0.48, respeetively, as 
one ehanges 62 from —0.20 to 2.00. On the other hand, when eorreetions and dimension-eight 
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0.2 
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0.0 0.5 1.0 1.5 2.0 
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FIG. 4. (Color Online) The ratios between quasi nucleon self-energies and the vacuum mass as a function 
of 62 - (a) shows the results without and (b) with as corrections. 
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0.00 
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0.05 
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-0.05' . -0.05 
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62 


62 


FIG. 5. (Color Online) Nuclear symmetry energy as a function of 62 without (a) and with (b) twist-4 spin-2 
four-quark matrix elements. The vertical axis is shown in units of GeV. For the detailed definition of the 


vector and scalar contributions, we refer the reader to IlSll . 


condensates are included, the corresponding values of and M^/Mpj are, respectively, 

modified to vary from 1.00 to 0.93 and from 0.58 to 0.52 for the same range of 62 - If 62 can be 
restrieted to the narrow range of Eq. (l26l) . this means that the dependenee of both and 

within this range elearly beeomes very small. Moreover, previous sum rule ealeulations 
based on the larger value of 62 given in Eq. (fT4l) [Il4n . have presumably overestimated the absolute 
changes of the above ratios in nuelear matter. 

Einally, let us eheck the modification of the symmetry energy as we change 62 as above. In 
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Fig. [51 we show the symmetry energy obtained with the full sum rule including correction and 
dimension-eight condensate without (a) and with (b) the four-quark twist-4 spin-2 contributions. 
As can be seen in the figure, the curves are essentially constant for both cases and hence show no 
dependence on 62 - This simply shows that, in contrast to the nucleonic parameters shown in Fig. 
m the contribution of 62 to the symmetry energy is negligible. 


IV. SUMMARY AND CONCLUSION 


We have in this work studied the matrix element ST{N\qDpiDvq\N) and have for the first 
time provided an estimate for its value that is based on experimental constraints. To do this, we 
partially had to rely on a number of assumptions on the relative contributions of u and d quarks and 
their respective antiparticles to 62 and, furthermore, on the behavior of e{x) at small and large x. 
For setting up these assumptions we followed the guidance of the bag model [|5|], the chiral quark 
soliton model [|9|] and the spectator model lilOll . which have been used to compute the contributions 
of the different quark flavors to e{x). 

The final result, given in Eq. (l26l) . contains both the uncertainties of the experimental results 
of [j^ and the systematic errors due to the assumptions derived from the different quark models. 
We should stress here that, even though uncertainties of Eq. (l26l) are large, the smallness of 62 
compared to the earlier estimates of Eqs.® and (IT4l) appears to be robust and does not depend on 
the details of our employed assumptions. It hence seems to be essentially impossible to reconcile 
the experimental data of [|4|] with the old estimates of 62 , which, therefore, should be discarded in 
future studies. 

To study the consequences of our findings, we have furthermore investigated the contribution 
of ST{N\qDpiDyq\N) to the OPE relevant for deep inelastic scattering and for the correlator of 
the vector meson current, that couples to the p, (0 and (j) meson states. Our calculations show that, 
with the new estimate of Eq. (l26l) . 62 is so small that its contribution to the OPE for all the above- 
mentioned cases turns out to be negligible. Eet us here especially mention the OPE corresponding 
to the (j) meson channel, for which the relevant operator is msST{N\qDp^Dyq\N). With the old 
value of Eq. ([T4l) . this would have become the dominant twist-4 term, with seizable consequences 
for the behavior of the ^ meson in nuclear matter, but our result shows that this is not the case. 
We furthermore studied the contribution of the present operator to the nucleon scalar and vector 
self-energy and the nuclear symmetry energy in nuclear matter as obtained from the nucleon sum 
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rule. From this, we found that the uncertainties of the self-energies due to e 2 are much reduced 
with our new estimate. For the nuclear symmetry energy, the effect of turned out to be minimal. 

Our new estimate of e 2 will also have consequences for the OPE of other channels not studied 
in this work, which will be relevant for investigations of the behavior of the respective hadrons in 
nuclear matter. One such case could, for instance, be the study of the D or B meson spectrum at 
finite density [I26ll27ll . As 62 , however, turns out to be small, it can generally be expected to have 
only a marginal effect, but this needs to be confirmed in a separate calculation for each channel of 
interest. The present authors have already started an investigation of the hyperon symmetry energy 
at finite density and found that it could indeed be non-negligible for this case 112811 . 
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